And Don't Forget The Black Holes 



H.A. Bethe^'* G.E. Brown ^ C.-H. Lee'^ 

^Floyd R. Newman Laboratory of Nuclear Studies, 
. . . Cornell University, Ithaca, NY 14853 

m ■ , 

[ Department of Physics and Astronomy, State University of New York, 

O ■ Stony Brook, NY 11794, USA 



O 



c 



Department of Physics, Pusan National University, Busan 609-735, Korea 
(E-mail: clee@pusan.ac.kr) 



Abstract 

> 

OO I The discovery of the highly relativistic neutron star (NS) binary (in which both 

■ NS's are pulsars) not only increases the estimated merging rate for the two NS's 

f-*) . by a large factor, but also adds the missing link in the double helium star model 

of binary NS evolution. This model gives ~ 20 times more gravitational merging of 
' low-mass black-hole (LMBH), NS binaries than binary NS's, whatever the rate for 

. the latter is. 

o 

• The recent discovery of Burgay et al. [1] of the double pulsar PSR J0737— 3039A 

ci • and PSR J0737— 3039B [2] is very interesting for many reasons. One that is 

not so obvious is that it involves just about the lowest possible mass main 
I sequence giants as progenitors, which have not been encountered in the evo- 

' lution of the other binary NS's, even though there are 20 to 30 times more of 

them, as we show below. 

The observational evidence is strong that NS binaries evolve from double he- 
lium stars, avoiding the common envelope evolution of the standard scenario 
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[3]. (A helium star results in a giant when the hydrogen envelope is lifted off 
- in binary evolution by being transferred to the less massive giant in the bi- 
nary.) In the standard scenario of binary NS formation after the more massive 
giant transfers its hydrogen envelope to the companion giant, the remaining 
hehum star burns and then explodes into a NS. In about half the cases the 
binary is not disrupted in the explosion. The NS waits until the remaining 
giant evolves (and expands) in red giant following its main sequence hydrogen 
burning. Once the envelope is close enough to the NS, the latter couples to it 
hydro dynamically through gravity. In the system in which the NS is at rest, 
it sees the envelope matter coming at it. Some of it is accreted onto the NS, 
although most flies by in the wake, being heated in the process, and is lost 
into space. The energy to expel the matter comes from the drop in potential 
energy as the orbit of the NS tightens. Formulas for the tightening and the 
amount of mass accreted by the NS were given by Bethe & Brown [4] . 

Chevalier [5] first estimated that in the common envelope evolution the NS 
would accrete sufficient matter to evolve into a black hole (BH). This was 
made quantitative by Bethe & Brown [4] who calculated that in a typical 
case, the NS would accrete ~ IMq, taking it into an ~ 2AMq low mass black 
hole (LMBH), similar to the BH we beheve resulted from SN1987A. 

Thus, if the NS had to go through common envelope evolution in a hydrogen 
envelope of ^ lOM© from the giant, it would accrete sufficient matter to go 
into a LMBH. Therefore, when the giant evolved into a helium star, which later 
exploded into a NS, a LMBH-NS binary would result provided the system was 
not broken up in the explosion. (About 50% of the time the system survives 
the explosion.) The above scenario was estimated [4] to take place 10 times 
more frequently than binary NS formation which required the two stars to 
burn helium at the same time, and, because of the greater mass of the BH, 
the mergings of binaries with LMBH to be twice as likely to be seen as those 
with only NS's. This is the origin of the factor 20 enhancement of gravitational 
mergers to be observed at LIGO, over the number from binary NS's alone. 

We review the known NS binaries and show that they are consistent with the 
two NS's in a given binary having very nearly the same mass, as would follow 
from the double helium star scenario. 

We list in Table 1 the 5 observed NS binaries with measured masses. The 
very nearly equal masses of pulsar and companion in 1534+12 and 2127+llC 
is remarkable. We show below that 1913+16 comes from a region of giant 
progenitors in which the masses could easily be as different as they are. 
The uncertainties in 1518+49 are great enough that the masses could well 
be equal. The main point in our letter is to show that the masses in the 
double pulsar J0737— 3039A and J0737— 3039B were probably very nearly 
the same before a common envelope evolution in which the first formed NS 
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Table 1 

5 observed neutron star binaries with measured masses. 



Object Mass (M©) Object Mass (Mq) Refs. 



J1518+4904 


i.OD_g 44 


J1518+4904 companion 


-■-■u^^-o.ii 


[6,7] 


B1534+12 


1 0009+0.0010 

l.DOOZ_Q QQ-^Q 


B 1534+ 12 companion 


1 Q4c;9+0.0010 
±.o^oz_o 0010 


[8] 


B1913+16 


i.44U8_0 0003 


B1913+16 companion 


1 0070+0.0003 
l.OO(O_0.0003 


[9] 


B2127+11C 


1 Q4Q+0.040 

1.041:y_Q Q4Q 


B2127+11C companion 


1 QfiQ+0040 
±.OUO_o 040 


[10] 


J0737-3039A 


1 007+0.005 
±.00/ _Q 


J0737-3039B 


i.ZOU_o 005 


[2] 



J0737— 3039A accreted matter from the evolving (expanding) hehum star pro- 
genitor of J0737— 3039B in the scenario of Dewi & van den Heuvel [11]. 



Since the hehum burning in the giant is an order of magnitude shorter than the 
hydrogen burning in time, for two giants to burn hehum at the same time they 
have to be very nearly equal in mass, within 4%. Thus the binary neutron star 
scenario is highly selective. If they do burn helium at the same time, they can 
go through common envelope evolution at this time. The hydrogen envelope 
of the slightly more massive giant expands in red giant and is transferred to 
the less massive giant, which in turn evolves in red giant. The time for these 
evolutions is so short that the helium stars are unable to accept the hydrogen 
which is lost into space [12] so that the final neutron stars are also close in 
mass. Each of the helium stars will explode, going into a NS. In about half 
of each explosion the binary is disrupted, but in about 1/4 of the cases it is 
preserved and a double neutron star results. 

Bethe & Pizzochero [13] used for SN 1987A a schematic but reahstic treatment 
of the radiative transfer problem which allowed them to follow the position in 
mass of the photosphere as a function of time. They showed that the obser- 
vations determine uniquely the kinetic energy of the envelope once the mass 
is known. Prom the envelope masses considered, the range of energies was 1 
to 1.4 X 10^^ ergs. Using the fact that the pressure following the supernova 
shock is radiation dominated, Bethe & Brown [14] showed from the known 
value of O.O75M0 of ^^Ni production, an upper limit on the gravitational mass 
of I.56M0 could be obtained. This is just the calculated Pe core mass for the 
known ISM© progenitor of 1987A. 

Calculations of Woosley are shown in Table 3 of Brown et al. [15] where the 
amounts of fallback material from distance of 3500 km and 4500 km are given. 
The fallback simply cancels the gravitational binding energy from the Pe core 
evolving into a NS. Since the density of matter is thin at the estimated 4000 
km bifuraction point, this cancellation is relatively insensitive to the precise 
point of bifurcation. Thus, we can use the Pe core mass as the mass of the NS. 

In Pig. 1, we show both recent and older calculations of Pe core masses. The 



3 



"T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — rn — I — I — I — I — r 



2.0 



„ 1.8 



03 

to 

CO 

1 1.6 
o 

CD 



1.4 



X^ 



• XX X 

_l I I I tX' I I I I I I I I I I I I I I I I I I I I I I I I L 

10 15 20 25 30 35 40 



M . , [M 1 
giant sun^ 



Fig. 1. Fe core masses for a grid of stellar masses. See text for explanation. 

filled circles and crosses correspond to core masses at the time of iron core 
implosion for a finely spaced grid of stellar masses [16]. The circles were cal- 
culated with the Woosley & Weaver code [17], whereas the crosses employ the 
vastly improved rates for electron capture and beta decay [18]. 

A rapid increase in Fe core masses occurs at Mgiant ~ 18Mq, just the mass of 
the progenitor of 1987A which we believe went into a LMBH. The Burrows & 
Woosley [19] evolution of 1913+16 took place in terms of Mgiant ~ 20Mq. (In 
this evolution the Fe core mass will be shghtly less than given in Fig. 1 because 
the mass transfer in binary evolution removes the hydrogen envelope, leaving a 
"naked" helium star. The wind loss from naked helium stars is great.) Now one 
can see that the difference in masses between, 19 and 20Mq giants is greater 
than O.IMq, so that the O.OSMq difference in masses between the pulsar and 
its companion in 1913+16 can easily be furnished out of the calculated Fe core 
masses. 

There is some tendency for the Fe core masses to decrease from Mgiant of 15 
to I6M0, where we would estimate 1534+12 and its companion to come from. 
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(The companion is slightly more massive than the pulsar.) 

The pulsar and companion in the double pulsar should come from Mgiant ^ 
IOM0, at the 1.25M0 of the present J0737-3039A. Following Dewi & van den 
Heuvel [11] , the first pulsar formed is increased in mass with mass transfer from 
the companion when the latter is a helium star. Low-mass helium stars, 2.3 
to 3.3Mq have to burn hotter than higher mass ones because of the greater 
energy loss from their surfaces. In reaction to this, they expand in a helium 
red giant phase. The pulsar now goes through a common envelope stage with 
the envelope of the helium star once the latter expands far enough to make 
contact. 

As noted, the NS accretes some of the matter and most is expelled in the 
wake. Thus, from the coefficient of dynamical friction Cd — 2 ln(6inax/&min) — 6 
[20], it is found [4] that 



We take Muca = 2.5Mq and MHc,f = I.SMq. Now the orbit in the double 
pulsar is decreased a factor of ~ 3.7^^^ from that of the binary NS. (The 
average of periods of 1913+16 and 1534+12 is 0.37 days, whereas that of the 
double pulsar is 0.1 day.) The 2/3 exponent is the Kepler relation between 
orbital period and radius, and we then find 



= 1.075 (2) 

-^pulsciTji 

which gives the correct 1.34M0 for J0737-3039B given Mpuisar,i = 1.25M0, 
the same as the unrecycled neutron star. The I.5M0 He star will burn into 
the unrecycled pulsar, with mass decreased ~ 15% by the general relativity 
binding energy. 



Therefore, in the double pulsar, the masses of the two giant progenitors must 
be very close, since the two NS's would have had nearly equal masses (as in 
1534+12), except for the mass transfer during the hehum red giant stage. 

The Salpeter mass distribution of giants gives the relative number propor- 
tional to (Mgiant) so for two equal mass giants, the distribution goes as 
( Mgiant )~^'^- Thus, the probability of two IOM0 giants is 26 times greater than 
that of two 2OM0 giants, as would be progenitors of 1913+16. Burgay et al. 
[1] estimated the ratio Noi^i/Nigi^ to peak at ~ 6, whereas Dewi & van den 
Heuvel [11] estimate that there are many of the PSR J0739-like systems not 
seen because of their weak signals compared with those from 1913+16, possi- 
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bly giving a factor ~ 30 in mergings. Our estimate based on the large number 
of low-mass giant progenitors confirms this. 

With ^ IOMq hydrogen envelopes and a proportional tightening of the orbit in 
the common envelope evolution in the standard scenario, the Bethe & Brown 
result of accretion of ~ IMq onto the first NS seems reasonable. After only 
part of this is accreted, the neutron star, like in 1987A, goes into a black hole. 

Why haven't we seen any LMBH-NS binaries? Van den Heuvel [21,22] has 
pointed out that NS's form with strong magnetic fields 10^^ to 5 x 10^^ gauss, 
and spin down in a time 

Tsd ~ 5 X 10^ years (3) 

and then disappear into the graveyard of NS's. (The pulsation mechanism 
requires a minimum voltage from the polar cap, which can be obtained from 
Bi2/P^ ~ 0.2 with Bi2 = B/IO^^G and P in seconds [22].) The relativistic 
binary PSR 1913+16 has a weaker field B ~ 2.5 x 10^° gauss and therefore 
emits less energy in magnetic dipole radiation. Van den Heuvel estimates its 
spin-down time as 10^ yrs. There is thus a premium in observational time for 
lower magnetic fields, in that the pulsars can be seen for longer times. 

Taam and van den Heuvel [23] found empirically that the magnetic field of a 
pulsar dropped roughly linearly with accreted mass. This accretion can take 
place from the companion in any stage of the evolution. A pulsar that has 
undergone accretion is said to have been "recycled" . 

Now as mentioned, in 1913+16 the pulsar magnetic field is ~ 2.5 x 10^*^ gauss 
and in 1534+12 it is ~ 10^° gauss. In J0737-3039A it is only 6.3 x 10^ gauss. 
These "recycled pulsars" will be observable for 100 times longer than a 
"fresh" (unrecycled) pulsar. 

The same holds for LMBH-NS binaries. The NS is certainly not recycled, so 
there is an about 1% chance of seeing one as the recycled pulsar in a binary 
NS. But we propose 10 times more of the LMBH's than binary NS's, of which 
we observe 5. Thus the total probabihty of seeing the LMBH binary should 
be about 50%. However, there may be additional reasons that the LMBH-NS 
binary is not observed. NS's in binaries in which the nature of the companion 
star is unknown are observed. We may have to wait for LIGO which will be 
able to measure "chirp" masses quite accurately. The chirp mass of a NS binary 
should concentrate near I.2M0, whereas the LMBH-NS systems should have 
a chirp mass of I.6M0, and there should be ~ 20 times more of the latter. 

In order to estimate the number of mergings per year for the LIGO we begin 
from the estimate of Kim et al. [24] and multiply by a factor of 26 to account for 
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the increase in this estimate by the large number of low-mass giant progenitors, 
neutron stars from two of these having been observed in the double pulsar. 
Additionally we multiply by the factor of 20 to include NS-LMBH binary 
mergings. This brings our estimate up to 0.5 to 3.6 mergings per year for the 
initial LIGO, and about 6000 times greater for the advanced LIGO. 
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